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The progress on the lattice computation of low moments of both the unpolarised and polarised nucleon structure 
functions is reviewed with particular emphasis on continuum and chiral extrapolations and comparison between 
quenched and unquenched fermions. 
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1. INTRODUCTION 

Deep Inelastic Scattering (DIS) experiments, 
such as eN — > eX or vN — > X form an im- 
portant basis for our knowledge of the structure 
of hadrons. In these processes the current probe 1 
(either a neutral current, j/Z , or charged cur- 
rent, W + /W~) with large space-like momentum 
— q 2 = Q 2 breaks-up the nucleon. The (inclusive) 
cross section is then determined by the structure 
functions Fx, F2 when summing over beam and 
target polarisations and, in addition, F3 when us- 
ing neutrino beams, and g±, gi when both the 
beam and target are suitably polarised. The 
structure functions are functions of the Bjorken 
variable x (0 < x < 1) and Q 2 . (Another class 
of structure functions - the transversity hi - can 
be measured, in principle, from Drell-Yan type 
processes or in certain semi-inclusive processes 
[2].) While the original pioneering discoveries 
were made over thirty years ago at SLAC, more 
recently experiments with polarised beams have 
been reported and the field remains very active. 
Recent experiments and proposals, [3,4], include 
HI and Zeus at DESY (unpolarised F2 at small 
x, [5] and F 3 , [6]), Hermes at DESY (polarised 
gi and g%, [7]), El 55 at SLAC (polarised gi, 32, 



*Plenary talk by R. Horsley at Lat02, Boston, U.S.A. 
1 A more complete set of structure functions available from 
DIS processes is given, for example, in [1]. 



[8]), Jefferson lab (structure functions in the res- 
onance region, [9,10]), COMPASS at CERN (po- 
larised gluon distribution, hi, A matrix elements, 
[11]), CCFR at Fermilab (unpolarised F 3 , [12]) 
and RHIC (spin physics, [13]). Recent results are 
given in the DIS conference series, [14]. 

A direct theoretical calculation of the structure 
functions seems not to be possible (but see [15- 
17]); however using the Wilson Operator Product 
Expansion (OPE) we may relate moments of the 
structure functions to matrix elements of certain 
operators in a twist or Taylor expansion in 1/Q 2 . 
Thus if we define 

q^H D ■■■iD q 
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then we have the Lorentz decompositions 2 

(N(p)\CF^-»™y - tr\N(p)) = 

2v n [p^---p^ -tr] 
(N(p,s)\0^ 5;{,7 ^-^ } -tr\N(p,s)) = 
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and the v„, a„, d„ and i„ can be related to mo- 
ments of the structure functions. For example we 
have for v n and F 2 
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dxx n - z F 2 {x 1 Q A ) = 

^4(f(^Q 2 ,n»«' )i5 M+o(i/Q 2 ) 



and similar relations hold between g\ and a n ; g 2 
and a linear combination of a n and d n ; hi and 
i n . Note that v n , a n (including the a n part of g 2 
- the so-called Wandzura-Wilczek contribution) 
and t n correspond to twist-2 operators and have a 
partonic model interpretation 3 ; d n is twist-3 how- 
ever, and does not have such an interpretation. 

Although the OPE gives v n from F\ (or F 2 ) for 
even n = 2, 4, . . .; v n from F 3 for odd n = 3, 5, . . .; 
a n from gi for n = 0, 2, . . .; a„, d n from 52 for n = 
2,4, . . ., other matrix elements can be extracted 
from semi-inclusive experiments, for example a\ 
by measuring tt^ in the final state, [18]. 

The sum in the previous equation runs over 
/ = u, d, s, g, . . .. We shall only consider f = u, 
d here and mainly the non-singlet, NS, or proton 
minus neutron (p — n) matrix elements when the 
f = s and g (gluon) terms cancel. These lat- 
ter terms are less significant for higher moments 
anyway as the integral is more weighted to x ~ 1 
when sea terms have less influence. The Wilson 
coefficients, E MS (1, g MS (Q)) are known perturba- 
tively (typically 2 — 3 loops). 

Present (numerically) investigated matrix ele- 
ments, [19,20], are v 2 = (x) (which may also be 
considered as a piece of the momentum sum rule 

E q (x) {Q) + (X) i9) =l)V3 = (X 2 ), V A EE (x 3 >, 

ao ~ Aq (with a connection to the quark spin 
component of the nucleon and also for Au — Ad = 
gA to the Bjorken sum rule), a\ ~ Aq^ 2 \ [21], 
a 2 - Aq (3 \ h - Sq, [22,23], t 2 - Sq (2 \ d x and 
d 2 . We shall mainly discuss here v n (n = 1, 2, 
3), ao, £o, ^i and d 2 . Earlier (lattice conference) 



3 Altcrnativc notations, based on the parton model are 



reviews include [24,25]. Since then emphasis has 
been placed first on results with 0(a) improved 
fcrmions, considerations of continuum and chiral 
limits, simulations with dynamical fcrmions and 
recently on the use of chiral fermions (which can 
ease the operator mixing problem) . Also possible 
higher twist contributions and ir, p and A ma- 
trix elements have been considered. We shall here 
briefly review progress in these fields. 

2. THE LATTICE APPROACH 

Matrix elements are evaluated on the lattice, 
[26], from ratios of (polarised or unpolariscd) 
three-point nucleon correlation functions to (un- 
polarised) two-point correlation functions, 



R a p{t,T;p) = 



(N(t;p)N(0;p)) 
as depicted in Fig. 1. Using transfer matrix meth- 
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Fi gure 1. The quark-line- connected diagram, left hand 
picture, and quark-line- disconnected diagram, right hand 
picture. The cross denotes the operator insertion 0{t). 

ods, it can be shown that R cx (N (p)\0\N (p)) 
provided that <C r <C t < \Nt (the lattice is of 
size Ng x Nt)- As the quark line disconnected di- 
agrams (RH figure of Fig. 1) are difficult to com- 
pute (for some reviews see [27,28]), it is again 
advantageous to look at non-singlet matrix ele- 
ments, such as v n; Ns = "i"' — = — v™. 
Finally most computations have been carried out 
in the quenched approximation, when the fermion 
determinant in the partition function is ignored. 
This is simply much cheaper in CPU time, but, 
as will be discussed later, unquenched results are 
beginning to appear. 

Although the Minkowski matrix elements dis- 
cussed in section 1 can be written in a Euclidean 
form in a straightforward way, the discretisation 
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onto a hypercubic lattice is not so restrictive as 
for the continuum and thus more representations 
appear, [29]. For example choosing the operators 

Oif b = q[j4D 4 - + l2 D 2 + l3 D 3 )}q 

both lead to a matrix element determining v 2 . 
The first representation requires a moving nu- 
cleon, while for the second a stationary nucleon is 
sufficient. An example for R for these bare matrix 
elements is shown in Fig. 2. Due to the increase 




Figure 2. The ratio R, normalised so that when the con- 
ditions 0< r<i < are met, the resulting plateau 
gives the bare matrix element. The picture shows the re- 
sults for quenched O (a) -improved fermions at (3 = 6.20 
and k = 0.1344 on a X N T = 24 3 X 48 lattice, t = 17 
and a typical fit range for the plateau is taken from r = 6 
to 11, shown by vertical dashed lines. The LH picture 
shows the 'diagonal' representation, i> 2 & with a station- 
ary nucleon, while the RH picture shows the 'off- diagonal' 
representation V2 a , which requires a moving nucleon with 
(lowest possible) momentum p = pi = (2iz/Ng,0,0). The 
empty symbols refer to u and d matrix elements, while the 
full symbols give the NS matrix elements. 

in noise in the signal, it is clearly advantageous 
to take a stationary nucleon - but this is unfor- 
tunately only possible for the lowest moments. 

The operators (or raw results) must be renor- 
malised. If using 0(a)-improved Wilson type 
fermions, one also wishes to improve the oper- 
ator. At present these additional operators are 
known for local (ie no D) and one-link (ie one 
D) operators [30], but not for higher- link oper- 
ators. Numerically when known these additional 



operators do not seem to be significant, [31]. Per- 
turbatively Z is known for all the local opera- 
tors, and for the one-link operators, [30]. For the 
higher-link operators only the unimproved results 
are presently known. For Ginsparg- Wilson (GW) 
fermions, the situation is simpler as the improve- 
ment coefficients are simple numbers, [32], while 
the renormalisation constants are given in [33]. 
Furthermore renormalisation constants for local 
operators for Domain Wall (DW) fermions have 
also been calculated in [34]. 

As most of the perturbative one-loop coef- 
ficients are known, then tadpole improvement 
(TI) of the renormalisation constants is possible. 
There are several variants, the one we shall use 
here is given in [35]. Here the renormalisation 
group invariant form of the renormalisation con- 
stant is directly computed. This can then be con- 
verted into the conventional MS-scheme. 

Non-perturbative renormalisation has also 
been attempted using both the Schrodinger 
Functional (SF) method, and the FJ-MOM- 
scheme. The SF method was developed mainly 
by the ALPHA collaboration, and presently in 
quenched QCD most of the improvement coef- 
ficients and renormalisation constants for 0(a)- 
improved Wilson fermions for the local opera- 
tors are known, [36,37], while for one-link op- 
erators (v 2a ) the renormalisation constant has 
been determined for both unimproved and 0(a)- 
improved fermions, [38]. An alternative ap- 
proach, RI-MOM, based on generalising the per- 
turbative procedure for the determination of the 
renormalisation constants has been applied to lo- 
cal, [39,40], and higher link operators, [40,35], for 
unimproved Wilson fermions. For 0{a) improved 
Wilson fermions local [41] and one-link operators 
[42] have been investigated. For DW fermions 
Z for local operators have been found in [43]. 
For unquenched fermions, very little is known at 
present, [44]. However as we later want to com- 
pare quenched and unquenched results, we shall 
use for consistency the TI Z (except for Z ao , [42] ) . 

Finally the operator mixing renormalisation 
structure is partially known. There is possible 
additional mixing with operators of the same di- 
mension for V3 and V4, [29]. Also mixing with 
lower dimensional operators occurs, in particu- 
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lar for d 2l [45], and d\. These are due to ad- 
ditional chiral non-invariant operators which oc- 
cur for Wilson fermions (but not for DW or GW 
fcrmions). This point will be discussed further in 
section 6.1. Also the formalism for the SF method 
has been developed for singlet operators mixing 
with gluon operators, [46]. 

3. CHIRAL AND CONTINUUM EX- 
TRAPOLATIONS 

The next step is to attempt a chiral extrapo- 
lation. An example for quenched 0(a)-improved 
Wilson fermions for v 2b is shown in Fig. 3. The 
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Figure 3. v%l s NS (2GeV) versus (rom ps ) 2 ~ rom q using 
(3 = 6.0, 6.2 and 6.4 for O (a) -improved quenched Wilson 
fermions. For orientation, the dash-dotted lines repre- 
sents (roughly) a strange pseudoscalar meson (m s quark 
mass, determined from rax ), an d the long-dashed line to 
the pion (ud quark mass). The chiral limit is given by the 
short dashed lines. Also shown is a linear extrapolation. 

points have been scaled to MS at \i — 2.0GeV 
(using r = 0.5fm ~ (400MeV) _1 ). Also shown 
is a linear extrapolation to the chiral limit, 

Vn-NS = a n (r m ps ) 2 + b n 

which seems to be adequate, but one should re- 
member that all the data points lie at the strange 
quark mass or higher. For w 3 and V4 similar ex- 
trapolations can be performed, but as they need 
a moving nuclcon yield much more noisy signals. 

Finally to obtain the phenomenological result, 
an extrapolation to the continuum limit must be 
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Figure 4. The continuum extrapolation for v^ s NS at 
a scale of 2GeV using the results from Fig. 3, filled cir- 
cles with the linear extrapolated result being given by the 
empty circle. Also shown is the LHPC+SESAM result, 
[20] using unimproved Wilson fermions, empty diamond. 
The stars are the MRS phenomenological values, [47]. 

performed. These extrapolations arc shown in 
Fig. 4 for V2, V3 and V4. One can see the degra- 
dation of the signal as one goes to the higher mo- 
ments, which makes the continuum extrapolation 
rather noisy. At least for V2b one can say that lat- 
tice effects appear to be small. Also shown is the 
result from [20] for unimproved Wilson fermions. 
Good agreement is seen, which again tends to sug- 
gest that lattice effects are small. The results of 
the extrapolation are also compared to the phe- 
nomenological MRS results, [47]. It is at present 
difficult to make any definite statement about the 
higher moments V3 and V4 except to say that due 
to the continuum extrapolation the ordering has 
been inverted (we would expect V3 > V4). The 
problem seems to lie in the continuum extrapola- 
tion and can only be cured with more /3-values. 
For V2 three [3 values seem to be sufficient; but 
the extrapolated result then seems to be about 
30% higher than the phenomenological value. 

One might be worried that one should perform 
the continuum extrapolation before the chiral ex- 
trapolation. The previous fits can be thought of 
as finding the best (3-(rom ps ) 2 plane to the data, 
so a variant procedure is to try a joint fit, [48], 

V2b:NS = a 2 (r m ps ) 2 + b 2 + c 2 (a/r ) 2 + d 2 ar m 2 ps 
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where the first two parameters represent the 
'chiral physics', the third parameter potential 
0(a 2 ) effects and the fourth parameter oc am q ~ 
arorrip S is to account for any residual quark mass 
effects. (With three (3 values, one actually re- 
duces the number of free parameters by one.) In 
Fig. 5 we show the results of this type of fit. The 
same continuum result is obtained. 
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Figure 5. v™ b s NS (2Ge V) °wiih lattice artifacts removed, 

ie V2b — C2(a/?"o) 2 — ^2<wo m ps showing the chiral extrap- 
olation for the quenched O (a) -improved Wilson results of 
Figs. 3 and 4- The same notation as in Fig. 3. 



4. UNQUENCHED RESULTS VERSUS 
QUENCHED RESULTS 

One possible explanation for the discrepancy 
between the lattice result for t> 2 and the phc- 
nomcnological result is the use of the quenched 
approximation. Indeed one might expect that 
due to the momentum sum rule the quenched re- 
sult is greater than the unquenched result (as the 
sea term part is suppressed in the quenched ap- 
proximation). While most of the data at present 
uses quenched fermions, some recent results us- 
ing unquenched fermions has appeared: from 
the LHPC and SESAM Collaboration, [20] (us- 
ing unimproved Wilson fermions with (3 = 5.5, 
5.6) and from the QCDSF and UKQCD Collab- 
oration (using 0(a)-improved Wilson fermions at 
(3 = 5.20, 5.25 and 5.29, [35,49]). (Both Collab- 
orations have three quark mass values at each (3 
value.) Again, as in the quenched lin- 
ear chiral extrapolation (at fixed (3) seems ad- 
equate. In Fig. 6 we plot the results against 
(a/ro) 2 . Are there quenching effects? Although 
the unquenched results are not as good as the 
quenched results, it seems that in this quark mass 
(> m s ) and a range quenching effects are small. 
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¥ • QCDSF quenched, O(a)-improved 
■ LHPC+SESAM quenched, unimproved 
) QCDSF+UKQCD unquenched,0(a)-improved 
□ LHPC+SESAM unquenched, unimproved 
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Fi gure 6. Quenched and unquenched results from the 
QCDSF+UKQCD and LHPC+SESAM Collaborations. 
Note that to determine a different scales have been used 
(either ro or rrajv being extrapolated to the chiral limit), 
and also slightly different renormalisation procedures. 

Further quantities that have been considered 
include the axial charge (ie Bjorkcn Sum rule) 

( dxgl~ n (x,Q 2 ) = \E gi:ao , NS 9A 
Jo 

(N(p,s)\q^ l5 q\N(p,s)) = 2s»Aq 

with Au^in) - Acft^in) = g A and the tensor 
charge 

/ dxh\ n {x,Q 2 ) — E™^ t NS tQ.% s 
Jo 

(N(p,s)\qi<J^ l5 q\N(p,s)) = —{s»p v - 3 v p^)5q 

with 5u~(fi) - Sd~(n) = t^sifi). In Figs. 7 
and 8 we show equivalent pictures to Fig. 6. 
Again little difference between the quenched and 
unquenched simulations is seen. 

5. TOWARDS SMALL QUARK MASSES 

The results shown previously have all been 
characterised by having data points at quark 
masses at or above the strange quark mass, and 
then a linear extrapolation in the quark mass to 
the chiral limit. There has been much recent work 
developing chiral perturbation theory, x _ PT, [50- 
53] which has shown the existence of a chiral log- 
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Figure 7. Quenched and unquenched results from the 
QCDSF+UKQCD and LHPC+SESAM Collaborations, 
for gA using the same notation as in Fig. 6. 



arithm of the form 



, In i 



tn:NS = 
V n ;NS = 
Cln;NS = 



T n (l-iC 4 m^ln(m^/^)) 

Vnil-CsmlMKs/ti)) 
A n (l- C 2 m 2 ps \n(mlJ^ 2 x )) 



where C n = {ng\ + l)/(47r/ 7r ) 2 for (full) QCD. 
For quenched QCD an expression for C3 in terms 
of F and D constants can be found in [54] . (Note 
that for quenched QCD for the nucleon there is 
no 'hairpin' contribution giving rise to the loga- 
rithm ~ \nm q for v n .) What is a suitable value 
for the chiral scale /i x ? Roughly for m ps > /i x , 
pion loops are suppressed and there is a linear 
variation in m q ie constituent quark behaviour, 
while for m ps < fi x we have non-linear behaviour. 
Often a value for yt x ~ lGcV is taken. We 
shall also use a comparison value of 500MeV here. 
From the above formulae, we see that X-PT al- 
ways decreases the value of the matrix element as 
m ps — > 0. Thus the lattice result should always 
be larger than the x-limit result. Also we would 
expect more effect for v n than for a n . (Recent 
work in [55] indicates however that when includ- 
ing A as well as the N then effectively bending 
only occurs for v n but not a n and t n .) 

As it is not so clear to which quark mass x~PT 
is valid, and as the results shown so far at m„ > 
m s yield a linear behaviour it is necessary to go 
to lower quark masses. In [35] this was started 
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Figure 8. Quenched and unquenched results from the 
QCDSF+UKQCD and LHPC+SESAM Collaborations, 
for (Su — Sd) MS (/i) using the same notation as in Fig. 6. 

using unimproved quenched Wilson fermions at 
(3 = 6.0 (as the problem of 'exceptional configu- 
rations' then seems to be less severe). The present 
status, [56], is shown in Fig. 9. All quark masses 




Figure 



Quenched results for unimproved Wilson 



fermions at (3 = 6.0 for v^ s NS (2GeV) . The empty circles 
are for heavier quark masses on a 24 3 X 32 lattice, O(200) 
configurations, while the filled circles for the lighter quark 
masses are on a32 3 x48 lattice, 0(100/50) configurations, 
[56]. Finite volume effects are checked at the second light- 
est mass. The dashed line is a linear fit to the data. From 
the MRS value (star), [47], the x-PT formula is applied, 
with Cs/rg ~ 0.28, 0.67 for quenched or full QCD respec- 
tively and for two values of the chiral scale /i x . 

have am ps Ns > 4. As (r 2m 7r ) 2 ~ 0.5, the light- 
est mass used in the simulation is ~ 2m u d (this 
corresponds to m ps /mv ~ 0.4). Little curvature 
in the numerical results is seen, but is still possi- 
ble, as we expect the coefficient C3 to be smaller 
than in the unquenched case, [54]. Also, as noted 
before, quenching might give a higher value for 
V2b than the phenomenological value anyway. 
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Fi gure 10. Quenched results for unimproved Wilson 
fermions at fi = 6.0 {or qa- Notation as in Fig. 9. 

In Fig. 10 we show results for gA- Note that \- 
PT goes in the wrong direction. It is less clear if 
there is a finite volume effect. Ref. [57] suggested 
that the charge is delocalised in the chiral (and 
infinite volume) limit, gA — > 2/3^. This was 
further discussed in [58], which showed that for 
finite volumes, there are no (large) volume effects. 

Nevertheless the questions of finite volume ef- 
fects and the range of applicability for Wilson 
fermions near the chiral limit remain, and re- 
cently there have also been results using DW 
fermions by the RBC Collaboration. These 
have much better chiral properties than Wilson 
fermions and so are more suitable for investigat- 
ing small quark masses. In [59], v^b is computed 
on a 16 3 x 32[xl6] lattice at a" 1 - 1.3GcV. Thus 
the lowest pion mass there, (am ps ) 2 <~ 0.1, corre- 
sponds to about (romp S ) 2 ~ 1 in Fig. 9. At this 
point some curvature in the signal is present. For 
gA, [60], finite volume effects are seen. Thus the 
general situation is not completely clear. 



6. OTHER TOPICS 

6.1. Non-perturbative mixing 

A further source of discrepancy between lattice 
results and phenomenological results can lie in the 
incorrect treatment of non-perturbative mixing of 
the lattice operators. An example is given by d%, 
which can be found from g<2, 



f 

Jo 

f 

Jo 



q—u,d 



dxx 2 gi(x, Q 2 



I X E ai"2 
q—u,d 



where ai and di are given in section 1 as nucleon 
matrix elements of certain operators C^/d 2 ~ 
QllnDDq. The a n operators have twist two, 
but d n corresponds to twist three and is thus 
of particular interest. A 'straightforward' lat- 
tice computation, [19,20], gave rather large val- 
ues for df. A recent experiment, [8], however 
indicated that this term was very small. This 
problem was traced in [45] to a mixing of the 
original operator O d q ^ with a lower-dimensional 

operator ~ quDq. This additional operator 
mixes oc 1/a and so its renormalisation constant 
must be determined non-pcrturbativcly. In [45] 
this was attempted using RI-MOM, and led to re- 
sults qualitatively consistent with the experimen- 
tal values. Note that this is only a problem when 
using Wilson-like fermions, as we would expect 
the operator to appear like ~ m q qcrDq and hence 
vanish in the chiral limit. Thus there should be no 
mixing if one uses GW or DW fermions. In [59] 
this was investigated for di using DW fermions 
and compared with unimproved Wilson fermions 
di results from [20]. The same phenomenon was 
seen: d\ using DW fermions gave a small value 
in the chiral limit, while the unimproved Wilson 
fcrmion results increased strongly as the quark 
mass was reduced. 

6.2. Higher Twist effects 

Potential higher twist effects are present in the 
moment of a structure function, see section 1. 
These 0(1/ Q 2 ) terms have four quark matrix ele- 
ments. A general problem is the non-perturbative 
mixing of these new dimension 6 operators with 
the previous dimension 4 operators. At present 
results are restricted to finding combinations of 
these higer twist operators which do not mix from 
flavour symmetry. In [61] the lowest moment of 
the pion structure function was considered, 

dxF2(x,Q )\Nachtmann = 

1.67(64) +0(a 2 s ) 

where the SUf{^) flavour symmetry group gives 
the combination F 2 /=2 = F? + + F?~ - 2Ff . For 
the nucleon the SUf{$) flavour symmetry group 
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must be considered, ie taking mass degenerate u, 
d and s quarks, [62], giving 

f dxF 2 {x 1 Q 2 )\ 2 ^= 1 tmann ^ 
Jo 

-0.0005(5)^^ + 0^) 

(To access this moment experimentally needs the 
measurement of structure functions of p, n, A, E 
and S baryons.) These results are for quenched 
unimproved Wilson fermions at j3 = 6.0, and are 
very small in comparison with the leading twist 
result. However these are rather exotic combi- 
nations of matrix elements and say little about 
individual contributions. Nevertheless this might 
hint that higher twist contributions are small. 

6.3. Pion, Rho and Lambda results 

Moments for pion and rho structure functions 
were computed in [63], for unimproved Wilson 
fermions. Using the SF method, v 2a was cal- 
culated for the pion, [64] for both unimproved 
and 0(a)-improved fermions, giving numbers in 
agreement with [63]. Finally there have been re- 
sults for moments of A structure functions, [65]. 
These are potentially useful as one can compare 
with nucleon spin structure and check violation 
of SUf(3) symmetry. First indications are that 
there is little flavour symmetry breaking. 

7. CONCLUSIONS AND FUTURE PER- 
SPECTIVES 

Clearly the computation of many matrix el- 
ements giving low moments of structure func- 
tions is possible. We would like to emphasise 
that a successful computation is a fundamental 
test of QCD - this is not a model computation. 
There are however many problems to overcome: 
finite volume effects, renormalisation and mix- 
ing, continuum and chiral extrapolations and un- 
quenching. At present although overall impres- 
sions are encouraging, still it is difficult to re- 
produce experimental/phenomenological results 
of (relatively) simple matrix elements (eg i>2, 9a)- 
Improvements are thus necessary in all areas. 
Nevertheless progress is being made: there are 
now considerations of both chiral and continuum 



extrapolations, some dynamical results are now 
available, there are attempts to understand lower 
quark mass both numerically (using both Wil- 
son fermions and DW fermions) and from x~PT. 
Clearly everything depends on the data and the 
quest for better results should continue. To leave 
the region where constituent quark masses give 
a reasonable description of the data, seems un- 
fortunately to require quark masses rather close 
to the ud mass. This will presumably also entail 
the use of unquenched chiral fermions (such as 
GW/DW). This will need much faster machines 
and is, perhaps, a cautionary tale for the deter- 
mination of other matrix elements. 
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